ON THE VACUUM POLARIZATION DENSITY CAUSED BY AN 

EXTERNAL FIELD 



CHRISTIAN HAINZL 

Abstract. We consider an external potential, — \ip, due to one or more nu- 
clei. Following the Dirac picture such a potential polarizes the vacuum. The 
polarization density, Pvaci as derived in physics literature, after a well known 
renormalization procedure, depends decisively on the strength of A. For small 
A, more precisely as long as the lowest eigenvalue, ei(A), of the corresponding 
Dirac operator stays in the gap of the essential spectrum, the integral over 
the density p^ ac vanishes. In other words the vacuum stays neutral. But as 
soon as ei(A) dives into the lower continuum the vacuum gets spontaneously 
charged with charge 2e. Global charge conservation implies that two positrons 
were emitted out of the vacuum, this is, a large enough external potential can 
produce electron-positron pairs. 

We give a rigorous proof of that phenomenon. 



1. Introduction 

In 1934 Dirac and Heisenberg realized that accepting the Dirac picture of elec- 
trons filling up the negative energy states, called vacuum, consequently implies that 
a charged nucleus thrown into the vacuum causes a redistribution of the Dirac sea, 
an effect denoted as vacuum polarization. Uehling and Serber in 1935 |271 I25|. 
long before standard renormalization procedure, demonstrated that such an indi- 
cated production of virtual electron-positron pairs give rise to a modification of the 
Coulomb potential and thus causes energy shifts of bound electrons. 

Concerning the traditional Lamb shift, known as the splitting of the 2sx/a- and 
2pi /2-state in hydrogen, this effect only accounts for about 2.5 percent. However the 
Uehling potential represents the dominating radiative correction in muonic atoms 
which emphasizes the importance of vacuum polarization (VP). Notice, whereas 
interaction with a photon field can be treated non-relativistically there is no non- 
relativistic equivalence for VP. It is a purely relativistic effect. 

Within the framework of QED, VP is treated by means of perturbation theory 
as developed by Dyson, Feynman, and Schwinger. 

Only recently Hainzl and Siedentop demonstrated in that the effective one- 
particle Hamiltonian obtained from VP can be handled non-perturbatively and 
gives rise to a self-adjoint operator. The effective potential we gain is in fact the 
same as the physicists obtain after mass and charge renormalization (neglecting 
photon terms) and use to calculate the hyperfine structure of bound states. We 
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refer to [201 Section 4] for a nice review concerning the influence of VP on the Lamb 
shift of heavy atoms. 

The main goal of the present paper is to study the vacuum polarization density 
caused by an external field, i.e., by one or more nuclei. As foreseen by physicists, 
e.g., |10llU). the behavior of the density turns out to depend on the lowest eigenvalue 
of the corresponding Dirac operator. As long as this eigenvalue stays isolated the 
integral over the density vanishes, that means the vacuum stays neutral. But as 
soon as that eigenvalue touches the lower continuum the vacuum gets spontaneously 
charged, i.e., an electron, more precisely two electrons due to degeneracy of the 
"ground state" , are trapped in the vacuum and two positrons are emitted. In other 
words large fields can produce electron-positron pairs. Such a situation can be 
realized by heavy ion collision. 

1.1. Model. The free Dirac operator is given by 

(1) D°:=a--W + (3 

i 

in which a, j3 denote the 4x4 Dirac matrices. The underlying Hilbert space is given 
by $j = L 2 (T) with T = M 3 x {1,2,3, 4}. We pick units in which the electron mass 
is equal to one. We regard the case of one, or more, smeared nuclei with density 
n e L 1 (M 3 ) n L°°(R 3 ), non negative, and assuming f R3 n = 1. 

We remark that it is an experimental fact that the nucleus cannot shrink to a 
point. In fact a point nucleus creates instability if one includes polarization effects, 
as shown in |111 Section 3.5]. 

The corresponding electric potential reads 

(2) ^ = \-\- 1 *n. 
and the operator to be studied is given by 

(3) D Xip := D° - X(p, 

where A > is a parameter and can be thought of as «Z, a the fine structure 
constant, e :— —yfa the charge of an electron, and — Ze the charge of the nucleus 
(nuclei). In the following we want to allow any value of A. 

Due to smearing out the Coulomb singularity the case of large values of A does 
not influence the behavior of the essential spectrum as well as the self-adjointness 
as it would be in the case of the Coulomb potential. The following Lemma is well 
known, e.g., Weidmann |281 Theorem 10.37]. 

Lemma 1. Lettp=\ - I" 1 *n, n € L^M 3 ) n L°°(IR 3 ), non negative. Then, VA > 0, 
jj\ip _ jjO _ ^ j s self-adjoint with domain H 1 ^) and the essential spectrum of 
D Xv is given by 

(4) a css ( J D A ^) = (-co,-l]U[l,oo). 

Throughout the paper we will denote the spectrum of D Xtp by a(D Xv ) and ej(A) 
as the corresponding eigenvalues. 

The following is well known: For fixed A there is an infinite number of eigenvalues 
which accumulate at 1 and each ej(A) depends continuously on A. For small values 
of A all eigenvalues stay in the gap (—1, 1) of the essential spectrum. However, for 
each i one finds a Xi such that for liniA^Ai e iW = ~ 1> i- e -> the eigenvalue ej(A) 
dives into the lower continuum. We do not at all discuss what happens to the 
eigenvalues after reaching the continuum. In fact one knows from [2] that below 
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— 1 there are no embedded eigenvalues. However, the behavior of the eigenvalues 
after reaching (— oo, — 1] won't play any role. Our theorems only depend on the 
number of eigenvalues, counting multiplicity, that vanish in the lower continuum. 
Namely, due to our assumption if > 0, all eigenvalues are monotonously decreasing 
(this is a consequence of [221 Theorem XII. 13] and the fact that each eigenvalue has 
non-positive derivative). That means they will not reappear after having reached 
— 1. The fact that each eigenvalue reaches —1 for a large enough parameter can be 
seen by e.g., a Theorem of Dolbeaut-Esteban-Sere 

We will see that whenever an eigenvalue dives into the "sea of occupied states" , 
i.e., (oo, —1], a specific number of e _ e + pairs are created depending on the degen- 
eracy of the dived eigenvalue. 

1.2. Vacuum polarization density. As already mentioned above, according to 
Dirac the vacuum consists of electrons occupying the negative energy states of the 
free Dirac operator. If one puts a nucleus into the vacuum, then the electrons 
rearrange and one ends up with virtual electron-positron pairs. In other words 
the vacuum gets polarized, see e.g., page 257], for a picture describing this 
phenomenon, and |1 lj for a "mathematical" derivation of the vacuum polarization 
density, which follows the idea of the early papers in QED 12H1 El E] • For a 
review about the old fashioned way of QED we refer to |18| . 
The operator describing this polarization effect is given by 

(5) Q x *> := P Xip - P°, 
where 

(6) P X J := X( -^-i]{D X n- 

Physically speaking we project onto the occupied states of the Dirac sea. 

Remark 1. Notice, in the case that the lowest eigenvalue of D Xv , ei(A), is strictly 
positive, our definition is equivalent to |111 Equation (12)], apart from a minus sign 
which is chosen to adapt to the definition in the physics literature. 

Usually the first idea to define a density via Q X{p would simply be taking the 
diagonal of the Kernel. Unfortunately, the operator Q Xv is not trace class. The 
question how to extract from Q Xv a physically meaningful density was first posed 
in the 30-ies by Dirac 0] and Heisenberg ^2J and in more recent literature this 
procedure is known as charge renormalization (see e.g., jSJE])- As in |11| we use 
Cauchy's formula to express the Q Xv in terms of the respective resolvents (Kato 
[T3| . Section VI, 5, Lemma 5.6) 



q\ v = p x v _ p o = _ dvj 



(7) ^ 2irJ_ oc \D°-j + ir) D x f-~f + ir\ 

._ pAtp pO 

where 

1 1 f°° 1 

(8) P^-=~-^Z dn 



7 ' 2 2vr7_ 0O D^-j + ir)' 

with — 1 < 7 < ej(A), e^A) being the lowest isolated eigenvalue of D Xv . Notice, that 
the second equality in J7J) is a consequence of the fact that D d_ is holomorphic with 
respect to z in the complex strip between (—1, 1) and D \l,_ z between (— 1, ej(X)). 
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We decompose Q Xv into 4 terms: 



(9) 

where 
(10) 

Qi ■= 



XQi + A 2 Q 2 



A 3 Q 3 



^ Q<1 ! 



i r i i 



n 1 r a 1 1 1 

<?2 := TT / dl ? no i ■ y no I ■ V no i • ' 



1 

2^ 



^ / dr l 



1 

2n 



(It) 



1 



1 



L>° + ir] 9 D° + iif D° + iif D° + irj' 



-9 



<P- 



D° - 7 + irj D° - 7 + if] D x v - 7 + ir\ D° - 7 + i?]^ D° - 7 + ny' 



The first three terms we consider by means of its Fourier representation. A simple 
variable transform irj — > irj + 7 does not change the Kernel of the operators Q\ 
to Q3 which is the reason why we suppressed the 7 in the denominator. The first 
term is treated in detail in |1 II Section 3.2]. There, by a well known renormaliza- 
tion procedure following Weisskopf [2H| and Pauli and Rose [21], we extracted the 
corresponding physical density 



(11) 



p\{x) := e\T~ x 



\k\< 



(*)> 



where (see ^] Equation (21)]) 



C(k)/k 2 



(12) 



2 jo 
1 



dx(l - x 2 ) log[l + k 2 {l - x 2 )/4] 



M 2. /" 4, V1 + 4/F+1 4 5 

(1 ~^W 1+ P l0g 7TTW r T + ^"3 



which was first explicitly written down by Uehling |27j and Serber |25j and later 
by Schwinger [22] and others (see also [H E3 OH]). Observe in [TT] Eq. (52)] that 
renormalization consists of subtracting an operator with infinite diagonal from Q±. 
From what rests one defines, in l|ll|) . the diagonal corresponding to Q±. This 
subtraction reflects the main difficulty concerning the proof of our main Theorem. 

The second and third term in l|10(l have a well defined integrable diagonal when 
using the Fourier representation. Additionally the density corresponding to Qi 
vanishes, either through integration over r\ or due to the fact that the Dirac matrices 
are traceless. Quite generally, if we expand tiQ-iQ^ into an infinite sum, each term 
with an even number of (p vanishes. 

The density corresponding to Q3 is given by 



(13) 



e\ 3 {2n)- 3 f dp f dqJ^^'^Q 



a(p,a;q,<r) 
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where Q3 denotes the Kernel of the Fourier representation 

(14) Qz(p,q) = ^-( dr] I dpx I dp 2 (D p + ir^y 1 o fi(p - Pl ) o (D P1 + i^y 1 

^J-oo JR 3 JR3 

o ip( Pl - p 2 ) o (D P2 + irj)~ l o 0(p 2 -q)o (D q + iriY 1 

with D r := a. ■ r + (3. 

The operator Q\ will be shown to be trace class in Lemma so we can define 
P4 quite general via the diagonal of X^tr^Q^, 

(15) p\(x) := e\h^Ql{x,x) 
Therefore the renormalized density reads 

(16) pt c (x) := f$(x) + p$(x) + pl{x). 

Before formulating our main theorem it is necessary to introduce the counting 
function d(\) which counts the number of eigenvalues which dived in the lower 
continuum for parameters smaller equal A. 
(17) 

d(X) :— {^eigenvalues, with multiplicity, that reached — 1 for parameters A < A} 
Theorem 1. Let n e L 1 (]R 3 ) n L°°(R 3 ) and non negative, ip = n * rj. Then 

(18) / p^ c (x)dx = ed(X). 



Theorem Q exactly reflects the picture which is presented by physicists, e.g., 
Greiner et al. piTU|: 

As long as the external potential, respectively A, is so weak that the lowest 
eigenvalue, ei(A), of D Xv is in the gap (—1,1) the vacuum stays neutral (and 
consists only of virtual electron-positron pairs). As soon as the lowest eigenvalue 
dives into the essential spectrum, (— 00, — 1], i.e., the sea of occupied states, the 
vacuum immediately gets charged with charge 2e (assuming that the ground state 
energy of D x,p is twice degenerate, due to the spin). This can be interpreted in the 
following way: when the unoccupied bound state dives in the sea of occupied states 
it traps two electrons which stay in the potential well of the nucleus (nuclei) . Due 
to Dirac's picture two "holes" emerge which are repelled and emitted as positrons 
out of the vacuum. Consequently we end up with real electron-positron, e~e + , pairs. 

This effect of spontaneously emitted positrons is verified in experiment by col- 
lision of heavy nuclei, which when approaching each other create an effective field 
strong enough to let the lowest eigenvalue dive into the continuum (see |22j ). 

Remark 2. In more recent physics literature, compare e.g., Equation (7.23)] or 
PUl Equation (230)], the VP-density is "formally" denoted as the diagonal of the 
operator 



e 



we 



(19) -trc^-P^L 

with := 1 - Since tr C 4[P« - Pi] = and -P^ + = P v - P« 

see that (|19H coincides with our initial operator etrc^P^ — P°_]. 

The proof of Theorem ^ will mainly be based on two ingredients: A work of 
Avron, Seiler, and Simon \I\ concerning the index of pairs of projectors (see also 
[7j) and arguments of Kato fTl) . 
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The proof of Theorem ^ will be given in Section |21 In Section we show that 
for tr[P V - p0]2m+l ) m > i ) a resu lt similar to JTSJ holds. 

2. Result on tr[P v - popn+i^ WITH m > j 

Recall that the vacuum polarization is in fact described by the operator Q Xip 
I ' — P°. Renormalization is inevitable, since that operator is not trace class. 
Nevertheless, due to Klaus and Scharf it is at least an Hilbert-Schmidt operator. 
Due to PP (in fact this follows already from Effros [7]) the traces of (Q Xv ) , 
m > 1, are equal. Therefore it is self-evident to ask for their behavior. 

Theorem 2. Let n e L 1 (R 3 ) n L°° (R 3 ) , non negative, tp = n* X. Then, Vm > 1, 

(20) tr[P^ - P ] 2 "^ 1 = d(X), 

where d(X) is defined as in 117JI . 



Proof. Notice that, since <p(k) = n(k)j§, 

(2D U gMtBtt < |tf| 

Jr3 1 + \k\ 



log(2 + |-|) 



H 2 (i + H) 



Take g = |, p = 4, then the second term on the right hand side is finite, as well as 
by Hausdorff- Young inequality 

(22) |||n| 2 || 4 = ||n|||<C7 8 3 /7 ||n||2 /7 <oo. 

Therefore the potential tp is regular in the sense of Klaus and Scharf > cf . |19l 
Equation (1.7)], namely the operator Q Xtp € & z ($j), i.e., Q Xv is an Hilbert-Schmidt 
operator. Consequently Q Xv e 6 m (S)) for any m > 2. 

To prove the Theorem we first look at the set of all A > such that the lowest 
eigenvalue, ei(A), corresponding to D Xip fulfills 

(23) ei(A) > -1. 

This is an open set so that we can always find a 7, with — 1 < 7 < ei(A) and 

(24) QXv = h L dr > - o^Tid = ^ - p " 

We are going to show that for m > 1 

(25) tr[P Xip - P°] 2m+1 = 

on the set {A|ei(A) > 7}. Since 7 can be chosen arbitrarily close to —1 this infers 
that 

(26) tr[P v - p0] 2m +! = 
on {A|ei(A) > -1}. 

To this aim we recall some results from Avron, Seiler, and Simon pQ (see also 
0) concerning the index of pairs of projections: 

Regard the family of orthogonal projections P Xv , A > 0. Since P Xv — £ 
62(8), P Proposition 3.2] implies that all pairs (P x T,PJfP) are Fredholm. Com- 
bining [Q Theorem 3.1] and P Theorem 4.1] we obtain that for m, I > 1 

(27) tr ^ ~ P ^ 2 " 1+1 = tr ^ ~ P ^ 21+1 = ind(P 7 V ' P ^ 

=dim(KerP^^ n RanP 7 Av ) - dim(KerP^ n RanP^) 
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is an integer. 

Remark 3. More generally, a pair (P, Q) of orthogonal projections is called Fred- 
holm, if the operator T = QP, as an operator from RanP — > RanQ, is Fredholm. 
The corresponding index ind(P, Q) is defined as 

(28) ind(P, Q) := indT = dim(KcrT) - dim(RanT) ± . 

Next we come back to the proof of i|25|) . Observe that on {A|ei(A) > 7} P Av 
is a continuous family with respect to the operator norm. Namely, using JHJ and 
expanding the resolvent we get 

r°° 1 1 

(29) ||P^-Pril<|A-^IIMI/ d V \\ 



f°° 1 



{51 + V 2 f 2 (S 2 2 +V 2 ) 1/2 ' 

with Si := min{7 + l,ei(/i) — 7} and 82 ■— min{7 + 1, ei(A) — 7}. Notice, that 
since ei(A) is continuous the integral in the right hand side of 129J1 can be bounded 
uniformly on a small enough closed neighborhood of each A in {A|ei(A) > 7}. 

Due to 1(23, \\P* V - P 7 W || < 1 implies \m\(P^,P^) = 0. Using Q] Theorem 
3.4(c)], 

(30) md(P^,P°) = ind(P A ^, i») + md(Pf , P 7 °), 

the continuity of P x,p immediately gives that ind(P Av ,P 7 ) = on the whole set 
{A|ei(A) > 7}. Together with (2ZJ we arrive at 

Summarizing, the argument given above was based on the fact that on the set 
{A|ei(A) > —1}, P Atp can be continuously deformed into PP. Throughout the rest 
of the paper we will repeat this argument several times. 

In the following we consider the case that an eigenvalue has dived into the lower 
continuum. We know that there are no eigenvalues below — 1. However, for nota- 
tional simplification we treat them as if they stay embedded. 

Fix now A such that ei(A) < —1 and 62(A) > — 1, and 7 with — 1 < 7 < 62(A). 
Additionally we choose a A' < A such that —1 < ei(A') < 7 and a 7' with —1 < 
7' < ei(A'). We know 

/ - \ 2m+l 

(31) tr (Q x *) = tr[P A ^ - pOfm+i = md (p^ p o^ 
Due to Theorem 3.4 (c)] 

(32) ind(P A ^ P 7 °) = md(P Av , P A '^) + ind(P A '^, P A >) + ind(P A >, J*,)- 

The first and third term in the right hand side in Q32JI vanish which can be seen 
by repeating the argument given above. Namely due to our choice of parameters 
P Xip can be continuously deformed into P A v . As well P , v can be continuously 
deformed into Py , which equals PP. 

Concerning the second term in the right hand side of l|32(l we note that by 
Cauchy's formula we obtain 

(33) P 7 Vv -P 7 > = P ei( v) 
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where P ei (\') is the projector on the eigenspace corresponding to the eigenvalue 
ei(A'). Consequently 

(34) ind(P 7 A >,P 7 A >)-tr[P ei(v) ]. 

By means of our definition l|17|l of d(\), obviously tr[P ei ( A /)] = d(X), whence 

/ t \ 2m+l 

(35) tr [Q x ^j = d(X), 



Repeating this argument whenever an eigenvalue dives into the lower continuum, 
(— oo, — 1], we arrive at the statement of the theorem. 

Notice, due to continuity in A the argument works no matter how many eigen- 
values "meet" at —1. □ 

3. Proof of Theorem [I] 

Summarizing the proof of Theorem [21 we exploited the fact that P_ ¥ ' build a 
continuous family of projectors on the non connected intervals 

(36) [0,A 1 )U(A 1 ,A 2 )U---U(A l ,A 4 +i)... 

where Xi denotes parameters where an eigenvalue reaches —1. As long as A, /z belong 
to a connected interval the index of the corresponding projection vanishes, 

(37) md{P^,P^) = 0, 

but if A moves to a different not connected interval the index jumps by an integer 
value. 

In order to prove Theorem \I\we first recall the definition of the density 

(38) pt c (x) = p$(x) + p$(x) + p$(x), 

the terms on the right hand side being defined in I|13ll . (|15|l . By means of our 
explicit choice of p* via Fourier transform Pi(k) — eXinh^r- and the fact that 
lim|fc|^ C ffi — we immediately obtain 

„A/ \ . -A, 



(39) / p\(x)dx = ft(0) = 0. 

Jr 3 

Therefore, our goal in the following will be to show that for all A 

(40) / pl(x)dx = 0, / pl{x)dx = ed(X). 
Jm 3 Jr 3 

It still remains to show that Q\ is trace class which works analogously to ^fl 
Lemma 3] . 

Lemma 2. 

(41) tr\Q X i \ = \\Ql\\i<C»\\y\\\, 

with an appropriate constant C M depending on // := min{7 + 1, e»(A) — 7}, &i{X) 
denoting the lowest isolated eigenvalue of D Xv . 
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Proof. Let ej(A) be the lowest isolated eigenvalue of D Xv> , then as usually, we choose 
a 7 with — 1 < 7 < ej(A). Using (|10J) we obtain (apart from a factor 

(42) ||Qi||i< 



< 



di) 



11111 

D° - 7 + irj^ D° - 7 + irj^ D x v - 7 + irj^ D° - 7 + irj^ D° - j + ir/ 
11111 



_D° — 7 + ir/^i) — 7 + ir/^ D° — 7 + irj^ D a — 7 + iij^ D a — 7 + it/ 

1 



-v- 



-<£>- 



( J D°-7 + i? ? )- 



with || (D° - 7 + iry) 
with 

(43) 



D*<p—~/+iri II — * 1 "IIV _dAvj_ 7 

with ||(£> +»„) I—I) < l + 



£)A V _ ^ _)_ ? ; r? 

|| < 1 + A||</3 DA *_ II which depends on Moreover, 



D°-7+i»7 1 
1 

-93 



1 



£)° — 7 + i?/ _D° — •-f + irj^ D° — 7 + z?/ 13° — 7 + iry ^ 13° — 7 + ?>/ 



— 



< 



rf?/ 



1 



1 



1 



1 



1 



3 



— V- 



(1 + M- 1 ) 5 



< 



1 



D° + it] 



1 



1 



<P 



D° + ir) D° + ir) 



1\5 



Applying an inequality of Simon |2fil Theorem 4.1], 

(44) \\mg(-W)\U < (27r)- 3 / 4 ||/|| 4 || 5 || 4 , 
to the factors in H43[l . gives 

llv^oTfelU * ^t^ti IMUIIWI- 1 2 + i + ^ 2 ll4 

(45) 1 x 1 

Putting all together and evaluating the integrals (cf. ^2 Lemma 3]) we arrive at 

(46) IIOjlli < ^imii, 

with an appropriate C M . □ 

In the following we will proceed analogously to the proof of Theorem [5] We will 
circumvent the problem that P_y — P° is not trace class by defining a family of 
trace class operators K E converging strongly to ip. We define K £ via its Fourier 
representation 

(47) 

with 



K e (p,q) :=f s (p)<p s (p-q)f s (q), 



(48) 



feip) '•= ~ |p|), <p e {x) := <p(x)x{l/e - \x\ 



X denoting the Heaviside step function. Obviously f e — > 1 and ip e — > y> pointwise 
when e — > 0. 

The family of operators 

(49) D XI<E := D° - XK £ 
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turn out to converge strongly to D Xip . For convenience wc define Q XK via an 
appropriate 7 chosen corresponding to D Xv , 

(50) Q AKe :=P 7 AKE -P°, 



leaving away the subscript 7, since it will not cause ambiguities. Furthermore 

DAK e _ p0 

■ 7 7 



P XK — P.P will be trace class so we can repeat the arguments given in the proof 



of Theorem |2 Since we already removed the "bad" part of P_ v — P°_ by charge 
renormalization (in 11 ), that is the part of Q\ which prevents P XK — P? from 
being trace class, it suffices to show that Q| (respectively Q\ ,e ) converge (in trace 
norm) to Q3 (respectively Q x ). Q% and Q^' £ are terms we obtain by expanding 

if5T7|. 

Notice, due to definition (fT3|) . J p^(x)dx — eA 3 J R3 trciQ 3 (p : p)dp. 
First we state a few useful properties of K e . 

Lemma 3. (a) For all e > 0, K £ is trace class and K £ > 0. Moreover, 
<t css (D xke ) = (-oo,-l]U[l,oo). 

(b) D XK ' -> D Xv strongly as e -> 0. 

(c) P XK ' - P° is trace class for all e > if J & a(D XI<e ). 

Proof, (a) The fact that K e is trace class is a direct consequence of Lemma [5] 
In Fourier representation we can decompose 

(51) K £ (p,q)=L* £ L E (p,q), 

with L e (p,p') = fs(p)h s {jp — p'), h e (p) = (27r) _3 / 4 y/^7(p). By our choice of f £ and 

(52) f f \L e (p,q)\ 2 dpdq<^ 

JR3 J R 3 

whence K e is trace class. Equation (|51(l immediately implies K s > 0. The com- 
pactness of K e yields a css (D XK ') = (—00, —1] U [1, 00) by Weyl's Theorem. 

(b) ForipeH 1 {T) 

(53) lim || [D XK ' - D Xv }t/j\\ = A lim || [K e - ^#|| = 0, 

~ * ~ 

since K 6 K e —> (p*ip in the sense of distributions, and these operators are bounded. 

(c) Let ef(A) be the lowest isolated eigenvalue of D XK . Then, with — 1 < 7 < 
4(A), 



(54) P^-P° = i-£ 



drj 



1 1 



D° — 7 + ir\ D XRe — 7 + in 
1 f°° 1 



2ttJ_ oc D°~j + in D XKs ~~f + in' 



Consequently 

1 f°° 1 

(55) \\P XK -POU^AII^H 1 — ^ — 

27r ^-oo ((l + 7) 2 + r; 2 ) 1/2 (8 2 +rf) ' 

with g ;= min{ 7 +l i e|(A)- 7 } using || D ^_ y+iv \\ < (P+ V 2 )- 1/2 and \\^^\\ < 

((i+^+^r 172 - □ 
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Let us fix an arbitrary A such that ei(A) > —1 and 7 with —1 < 7 < ei(A). 
Since D XK " -> D x f strongly, a Theorem of Kato [H VIII-5, Theorem 5.1] tells us 
that a{D XK ) is asymptotically concentrated in any open set containing a(D Xv ). 
Thus we can find a S small enough such that 7 < ei (A) — 6 and a corresponding £0 
such that for all e < Eq, cr(D XKe ) is concentrated in a <5-nei ghborhood of cr(L> Aip ), 
in particular ef(A) > ei(A) — 5 for each eigenvalue ef(A) of D XK \ 

Thus we are able to guarantee that P XK can be continuously deformed into P°. 
Therefore we can argue analogously to the proof of Theorem [S] combined with the 
trace class property of P XK — P° to obtain 

(56) tr[F 7 AKE - P°] = ind(F 7 AKE , P°) = 0. 

Expanding the resolvent this implies 

(57) 

r°° 1 1 f°° 1 1 1 

= Atr / dn-—-K £ -—-+\Hr / d„——K — —A 



D° + irj D° + irj J_ 00 D Q + irj D° + irj D° + irj 
A 3 tr / dn— K e — K £ — K e 



D° + ir) D° + irj D° + irj D° + irj 
A 4 tr / drj— jr—^ K £ ——^ K £ x 

1 -K^ K- 1 



D XKE — 7 + irj D° — 7 + irj D° — 7 + irj 

:= A tr Q\ + A 2 tr Q\ + A 3 tr Q% + A 4 tr Q A ' £ . 

Observe, this holds in particular in a small neighborhood of 0. Thus, since the 
fourth term on the right hand side is of order 0(A 4 ), each term in 1(5 7j) vanishes 
separately. In particular these terms we are interested in: 

(58) trQl = 0, trQ A ' £ =0, 

the latter one on the set {A|ei(A) > —1}. 

Assume for a moment we have already shown 

(59) eA 3 limtrQ! = / p X {x)dx, eA 4 lim tr Q A ' e = / p$(x)dx, 

then by 158(1 obviously L 3 p x (x)dx = and f R3 p x (x)dx — whence Theoremnon 
{A| ei (A)>-l}. 

In order to prove (|59|l we formulate an auxiliary Lemma: 

Lemma 4. (a) There exists a non negative function g £ L 4 (R 3 ), such that 

(60) |trc«£i(p,p)| <g(p) 
uniformly in e. 

(b) Let ei(X) be the lowest isolated eigenvalue of D Xlf> . Fix 7 and 5 with —1 + 5 < 
7 < 6j(A) — (5. Furthermore fix eo such that for all e < eo, <r(D XKe ) is in a 
5 -neighborhood of a {D x,p ). Then 

(61) IIQ4' £ lli<nM| 4 
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uniformly in e < Sq, where C M is an appropriate constant depending on )i 
min{7 + 1 — 6, e«(A) — 6 — 7}, and 



(62) 



lim trQ 4 ' e = tr t^ 4 . 



Proof, (a) We will proceed similarly to [111 Lemma 4]. For completeness we will 
repeat some parts of the proof. The "eigenfunctions" of the free Dirac operator in 
momentum space are 



(63) 



i T (p) ■= < 



a ■ pe T 
(l-E(p))e Ti 
a ■ pe T 
\-(l + E(p))e T 



r=l,2, 
t = 3,4 



with e T := (1,0)* for r = 1,3 and e T := (0,1)* for r = 2,4 and 

(64) N+(p) = y/2E(p)(E(p) - 1), N-{p) = j2E(p){E(p) + \). 

The indices 1 and 2 refer to positive "eigenvalue" E(p) and the indices 3 and 4 to 
negative —E(p). Using Plancherel's theorem we get 

4 ^4 

(65) tr C AQl(p,p) = V (u TQ {p)\Ql\u To {p)) = — — j / dpi dp 2 V x 

X (li To (p) I X £ I U T1 (pi )) (U T1 (pi ) I AT £ I U T2 (p 2 )) (li r2 (>2 ) I ^ £ I Ur (p)> 

1 



drj- 



{ia To E{p) - n^ia^Eip!) - r))(ia T . 2 E(p 2 ) - r))(ia To E(p) - rj) ' 

with a T = 1 for r = 1,2 and a T = —1 for r = 3,4. The integral over rj is seen 
to vanish by Cauchy's theorem, if all four a T - have the same sign. In fact we only 
treat one case. The others then work analogously. 
Set 



(66) 

Using f e < 1 the first factor in (|65|l can be estimated by 
(67) 



T 2 =3,4 

a-pa- Pl + (l-E(p))(l-E( Pl )) 



t = 1,2 ri=l,2 
< \<f>(p - Pl)<p(pi -p 2 )<f(p 2 -P)\ 



tiy 



N_(p 2 yN + {p) 2 N + (piY 
x [a ■ p x u ■ p 2 + (1 - E( Pl ))(l + E{p 2 ))} [a ■ p 2 a ■ p + (1 + E(p 2 ))(l - E(p))] 

\p-P2- (E( P2 ) -!)(! + E(p))\ + \ P A P2 \ 



< C\tf(p - pi)tp(p! -p 2 )lf>(p 2 ~p)\ 

(c being a generic constant.) Since 
1 f°° 

m T J_Jn 



N-(p 2 )N + (p) 



1 



(iE(p) - rf){iE{pi) - rj){-iE(jp 2 ) - f])(iE(p) - rj) 



2{E{p)+E{p 1 )fE{p) 
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our term of interest l|65(l is bounded by a constant times 



(69) 



dpi 



dp 2 \0(p - Pi)<p(pi -p 2 )(p(p 2 ~p)\x 

„ \p-P2- (E(p 2 ) -!)(! + E(p))\ + \ P A P2 \ 



2N^( P2 )N + (p)(E(p) + E{p x )fE{p) 
Substituting p 2 — > p 2 + P, Pi — » Pi + P2 + P we get 



(70) \m\< [ dpi [ dp 2 \0(p 1 +p 2 )0{pi)0(j) 2 )\'. 



\p ■ (P2 +p)- (E(p 2 +p) -!){! + E(p))\ + \ P A (p 2 + p)\ 
2N_(p 2 + p)N + (p)(E(p) + E( Pl +p 2 +p)) 2 E(p) 



Since 



\p ■ (P2+P) - (E(p 2 +p) - 1)(1 + E(p))\ + \ P Ap 2 \< A\p\\p 2 \ 
we obtain as an upper bound the function 

(71) g(p) ' ■ ' ' ' 



dpi / dp 2 \<p(pi + p 2 )<p(pi)<p(p 2 )\\p 2 
Jn 3 

which is obviously in L 1 (R 3 ), whence (a) is proven 



N_(p 2 +p)E{pY' 
(b) Analogously to (|42|l and (|43l) we get (apart from a constant) 



-( 



dr/x 



-K 



-K 



K 



-if 6 



D a - -f + ir] D° --f + ir) D XK " - >y + ir) D° - 7 + ir] D Q - 7 + ir] 



< 



dt) 



1 



-K £ 



1 



1 



1 



-K e 



1 



D° + ir] D° + ir] D° +ir) D° + ir\ D° + ir] 

1 



x 1 + 



D 



7 



The first term in the third line is trace class, so we can evaluate it in Fourier 
representation. Since f s < 1 we are in the situation of Lemma |21 and end up with 



(73) 



which implies l|61(l since \\K £ \\ is uniformly bounded. 
In order to prove H62J) it suffices to show 



\Q x 5 ' e - Q x 5 \\i 



0. 



since tr Q4 = 0, with 
1 



(74) Q 4 



2tt 



dr\ 



1 



D° + irj^ D° + irf D° + irf D° + irf D° + irj 



—V- 



-tp- 



In fact, due to Lemma |21 Qi is traceclass and using the Fourier representation one 
easily sees that its trace vanishes. Indeed each operator Q 2n with an even number 
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of potentials has vanishing trace, which is well known as Furry's Theorem. For 
convenience we denote 

i r°° 



(75) 



1 f°° 

Q^ = — drjN^ M v M n M v M v , 

27TJ-00 



where 2V2 = -p if 6 nu 1 , . , Ml = K e n0 1 ,. , 7V„ = n> „ 1 , . </? n o * 

and = ^ ^0^- ■ A straightforward calculation gives 



\N*M*M*Me\\ i/3 



1 f°° 

(76) ||Q^ 8 -Q^|| 1 <bu P ||M=-M 1j |U t - / ^ 

+ \\N*M*M*Mj A/3 + HA^M^IU/a + \\N^M n M v M n \U /3 

1 f 00 

+ sup \\m - AT,, a,— / A^M^M, i. 

Notice, since trace ideals fulfill 6 p C ©i, for p > 1, we can easily estimate both 
integrals in (|76fl in analogy to (|61fl and l|41|) . 
Obviously 



(77) ||M« - M„|| 4 < 



1 



vV+T 



dpdqf e {p, q) 



1/4 



Recalling the definition of K e in 1(47(1 . we see that f e (p,q) -^e^o pointwise, as 
well as 



(P,g) < TXT [101 * I^0 2 (P- «)^-T e l1 ( r3 x r3 ) 
p z + 1 5+1 



which implies, by the diminated convergence theorem, that sup^ ||Mf — M^||4 -^ e ^o 
0. 

Notice that 

( 78 ) n>yt 1 -. i .-_ -> FTx^ T , ~ strongly 



and 
(79) 



1 1 



K 



D° — 7 + ir] D a — 7 + if] 

both uniformly in r\. Together with the fact that if D iy_^ +iri is compact (even in 
64), i.e., it can be approximated in norm by a finite rank operator, we conclude 

sup||A^- Ag|oo -^0 
i) 

which yields 1)61(1 . □ 
Now we are ready to prove (jl)9"|> . Obviously, due to our definition il4"?|) and (|4*5|) , 

(80) tr C 4<3!(p,p) ^tr C 4Q 3 (p )P ) 

pointwise as e — > 0. By means of Lemma 01 (a) and the dominated convergence 
theorem 



in 64, 



(81) 



£^0 



eA 3 lim / tr C 4Q^(p,p)dp = eA 3 / tr C 4.Q 3 (p,p)dp 



R3 



R3 



p^(a;)da 
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By means of Lemma 01 (b) we obtain 

(82) eA 4 limtrQ,* 6 = eA 4 trQ^ = / p$(x)dx, 
whence 

Fix again A such that ei(A) < —1 and 62(A) > — 1. (As in the previous section 
we use the notation ei(A) for convenience. The argument works whatever happens 
to the eigenvalue after reaching the lower continuum. Let us remark, that for 
our results only the features of the eigenvalues before "diving" play a role.) We 
can find a 6 > small enough such that the following holds: We can choose a 
7 with —1 + 6 < 7 < 62(A) — S. Additionally we choose a A' < A such that 
— 1 + S < ei(A') < 7 — S and a 7' with — 1 + S < 7' < ei(A') — S. Moreover we find, 
due to Kato [141 VII-5, Theorem 5.1], an Eq such that for all e < Eq, er(P/ AK;e ) is in 
a ^-neighborhood of cr(D Xv ) as well as a(D x K ) in a ^-neighborhood of a(D x v ). 
We can write 

(83) tr[P 7 * Ke - P°] = tr[P 7 XKS - P^\ + - J$' K '] + tr[P^ - P°]. 

By our choice of parameters P 7 K ^ can be continuously deformed into P x K ^ , as well 
as Py K into Py which equals P®. Consequently the first and the third term on 
the right hand side of 1)83(1 vanish. Due to Cauchy's formula 

(84) tr[P 7 VKe -P^ K >tr[P ei(A0 ], 

due to the fact that by our choice of parameters the eigenspace of the set {e^(A')|7' < 
Bj(X') < 7} has the same dimension as the eigenspace of ei(A'). Recall P ei (A') de- 
notes the projector on the eigenspace corresponding to ei(A'). By definition l|17|) 



(85) tr[P ei(A0 ] = d(A). 
Whence 

(86) tr[P* K ° - P 7 °] = d(A). 

Expanding the left hand side as in l(57|l and using the fact that we already know 
that the first three terms vanish we see 

(87) A 4 trQ-p = d(X). 

By means of Lemma 01 (b) lim e ^o tr Q\' e = tr Q\, so we infer 

p x {x)dx = ed(X). 

Repeating this argument whenever an e^(A) dives into (—00, —1] we arrive at the 
Theorem. 



Appendix A. Criterion for a Kernel to be trace class 

It is well known that given an integral operator via a Kernel K(x,y), the fact 
that J R „ dxK{x,x) < 00 does not at all guarantee that K is trace class. 

For a specific class of Kernels we give a sufficient condition for the corresponding 
operator to be trace class. 
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Lemma 5. Let 

(89) K(x,y) = f 1 (x)g(x-y)f 2 (y), 
with /i,/ 2 € L 2 (R") and g € L 1 (M n ). Then K is trace class. 
Proof. We can write 

(90) g(k) = |g(fc)| 1 /2|5 (fc) |i/2 sgn( ^ (fc)) . 
Define 

(91) h x {x) := {2^)-^T- l [\g\ 1 ' 2 ] (x), h 2 (x) := ^[Isl^sgn^)]^), 
such that 

(92) g = T- 1 [g] = (2^- 3 / 2 ^- 1 [|. 9 r /2 ]*^- 1 [lffr /2 sgn(g)] = /ii*/ i2 . 
Therefore 

(93) K{x,y) = dzL 1 (x,z)L 2 (z,y), 



with 

(94) L 1 (x,z) = f 1 (x)h 1 (x-z) L 2 (z,y) = h 2 (z-y)f 2 (y). 
Observe 

(95) J dxdz\V{x, z)\ 2 = J J dxdz\f 3 {x)\ 2 \h 3 {z)\ 2 = \\f 3 \\ 2 2 \\gh < oo 
for j = 1,2, which implies the Lemma. 

□ 
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